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Magnetic insulators have proved to be fertile ground for studying new types of 

o 

quantum many body states, and I survey recent experimental and theoretical exam- 
pies. The insights and methods transfer also to novel superconducting and metallic 
states. Of particular interest are critical quantum states, sometimes found at quan- 
tum phase transitions, which have gapless excitations with no particle- or wave-like 
interpretation, and control a significant portion of the finite temperature phase dia- 
gram. Remarkably, their theory is connected to holographic descriptions of Hawking 
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radiation from black holes. 



I. INTRODUCTION 

o 

Quantum mechanics introduced several revolutionary ideas into physics, but perhaps the 
most surprising to our classical intuition was that of linear superposition. This is the idea 

> 

that any linear combination of two allowed quantum states of a system is also an allowed 
state. For single particles, linear superposition allows us to explain experiments like double- 
t-h slit interference, in which electrons can choose one of two possible slits in an obstacle, 

but instead choose a superposition of the two paths, leading to an interference pattern 

O 

in their probability to traverse the obstacle. However the concept of linear superposition 
had unexpected consequences when extended to many particles. In a well-known paper, 
Einstein, Podolsky, and Rosen pQ (EPR) pointed out that even with just two electrons, 
linear superposition implied non-local and non-intuitive features of quantum theory. They 
imagined a state in which the spin of one electron was up and the second down, which 
we write as |f|). Now let us take a linear superposition of this state with one with the 
opposite orientations to both electron spins, leading to the state (|tl) - |lt))A/2 Such 
a 'valence bond' state is similar to that found in many situations in chemistry and solid 
state physics e.g. the state of the two electrons in a hydrogen molecule. EPR imagined a 
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'thought experiment' in which the two electrons are moved very far apart from each other, 
while the quantum states of their spins remains undisturbed. Then we have the remarkable 
prediction that if measure the spin of one of the electrons, and find to be ||) (say), the 
perfect anticorrelation of the spins in the valence bond implies that the spin of the other 
far-removed electron is immediately fixed to be ||). EPR found this sort of 'action-at-a- 
distance' effect unpalatable. Today, we have many examples of the EPR experiments in the 
laboratory, and the peculiar non-local features of quantum mechanics have been strikingly 
confirmed. In recent years, the term 'entanglement' has frequently been used to describe 
the non-local quantum correlations inherent in the valence bond state of two spins. 

The purpose of this article is to describe the consequences of the linear superposition 
principle of quantum mechanics when extended to a very large number of electrons (and 
other particles). We are interested in how many electrons entangle with each other in 
the low energy states of Hamiltonians similar to those found in certain transition metal 
compounds, and related 'correlated electron' systems. While there is a precise definition of 
entanglement for two electrons which captures the essence of the EPR effect, much remains 
to be understood in how to define an analogous measure of 'quantum non-locality' for large 
numbers of electrons. We shall sidestep this issue by restating the question: how can 
we tell when two quantum states have genuinely distinct superpositions of the electron 
configurations ? Further, can we classify and list the distinct classes of states likely to be 
found in correlated electron materials ? 

A useful strategy is to follow the low energy states of a Hamiltonian which varies as 
a function of a parameter g. We examine whether the ground states breaks a symmetry 
of the Hamiltonian. We also examine the quantum numbers of the excitations needed to 
describe the low energy Hilbert space. If one or more of such 'discrete' characteristics is 
distinct between large and small g, then we assert that these limiting states realize distinct 
quantum phases. The distinction between the phases implies, in particular, that they cannot 
be smoothly connected as g is varied. Consequently there must be at least one special value 
of g, which we denote g c , at which there is a non-analytic change in characteristics of the 
ground state e.g. the ground state energy. Often some of the most interesting quantum 
states, with non-trivial and non-local types of entanglement, are found precisely at the 
quantum critical point g = g c , where the system is delicately balanced between the states 
of either side of g c . This happens when there is a correlation length which becomes large 
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near g c , as is the case close to a "second-order" phase transition (near strong "first-order" 
transitions the ground state simply jumps discontinuously between simple limiting states 
on either side of g c ). Thus a description of the Hamiltonian as a function of g exposes two 
quantum phases and their connection across a quantum phase transition, and we present 
several examples of such phenomena below. 

Experimentally, it is almost never possible to examine the ground state of a large quantum 
system in isolation. There are small random perturbations from the environment which 
eventually bring the system into equilibrium at a temperature T > 0. So it is crucial to 
examine the fingerprints of the ground state quantum phases and quantum critical points 
on the non-zero temperature phase diagram. We will turn to this issue in Section |Vj where 
we will argue that such fingerprints can often be significant and easily detectable. Indeed, 
the quantum critical point can serve as the best point of departure for understanding the 
entire phase diagram in the g, T plane. 

We will begin in Section [TT] by describing the rich variety of quantum phases that appear 
in two-dimensional insulators, whose primary excitations are S = 1/2 electronic spins on the 
sites of a lattice. Interest in such two-dimensional spin systems was initially stimulated by 
the discovery of high temperature superconductivity in the cuprate compounds. However, 
since then, such quantum spin systems have found experimental applications in a wide 
variety of systems, as we will describe below. The critical quantum phases and points 
of such materials will be discussed in Section IHIl Section [TV] will extend our discussion 
to include charge fluctuations, and describe some recent proposals for novel metallic and 
superconducting states, and associated quantum critical points. 

Although we will not discuss these issues here, related ideas apply also in other spatial 
dimensions. Similarly quantum states appear in three dimensions, although the quantum 
critical points often have a distinct character. A fairly complete understanding of one di- 
mensional correlated states has been achieved, building upon 'bosonization' methods. 

II. PHASES OF INSULATING QUANTUM MAGNETS 

High temperature superconductivity appears when insulators like La2Cu04 are doped 
with mobile holes or electrons. La2Cu04 is an antiferromagnet, in which the magnetic 
degrees of freedom are S — 1/2 unpaired electrons, one on each Cu atom. The Cu atoms 



4 



reside on the vertices of square lattices, which are layered atop each other in the three- 
dimensional crystal. The couplings between the layers are negligible (although this small 
coupling is important in obtaining a non-zero magnetic ordering temperature), and so we 
need only consider the 5 = 1/2 spins S», residing on the sites of a square lattice, i. The 
spins are coupled to each other via a superexchange interaction, and so we will consider here 
the Hamiltonian 

H = ^2 J ijSi ■ S i (!) 

(ij) 

Here (ij) represents nearest-neighbor pairs, and > is the exchange interaction. The 
antiferromagnetism is a consequence of the positive value of which prefers an anti-parallel 



alignment of the spins. We will initially consider, in Section |II A| the case of i on the sites 
of a square lattice, and all Jij = J equal, but will generalize in the subsequent sections to 
other lattices and additional interactions between the spins. 

Our strategy here will be to begin by guessing possible phases of H , and its perturbations, 
in some simple limiting cases. We will then describe the low energy excitations of these 
states by quantum field theories. By extending the quantum field theories to other parameter 
regimes, we will see relationships between the phases, discover new phases and critical points, 
and generally obtain a global perspective on the phase diagram. 



A. Neel ordered states 



With all = J, there is no parameter, g, which can tune the ground state of H ; J sets 
the overall energy scale, but does not modify the wavefunction of any of the eigenstates. So 
there is only a single quantum phase to consider. 

For the square lattice, there is convincing numerical and experimental evidence [2] that 
the ground state of H has Neel order, as illustrated in Fig. [T£i. This state spontaneously 
breaks the spin rotation symmetry of Ho, and is characterized by the expectation value 

(Si) = (-1)^ (2) 

Here (— 1) J represents the opposite orientations of the spins on the two sublattices as shown 
in Fig [pi. The vector <3? represents the orientation and magnitude of the Neel order. We 
have < 1/2, implying that the ground state wavevefunction is not simply the classical 
state sketched in Fig. [l^, but has quantum fluctuations about it. These fluctuations will 
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FIG. 1: Neel ground state of the 5 = 1/2 antiferromagnet Hq will all Jjy = J on the (a) square 
and (b) triangular lattices. The spin polarization is (a) collinear and (b) coplanar. 



entangle the spins with each other, but qualitatively the essential character of the state is 
captured by the pattern in Fig. [pi. Note that this pattern implies a long-range, and classical, 
correlation between the spins, but not a significant amount of entanglement because there 
are no EPR effects between any pair of well separated spins. 

Having described the ground state, we now turn to a consideration of the excitations. 
This is most conveniently done by writing down the form of the Feynman path integral for 
the trajectories of all the spins in imaginary time, r. After taking the long-wavelength limit 
to a continuous two-dimensional space, this path integral defines a quantum field theory in 
2+1 dimensional spacetime with co-ordinates (r, r). The quantum field theory is for a field 
<fr(r, r), which is the value of the Neel order when averaged over the square lattice spins, Sj 
located on sites within some averaging neighborhood of r. A quick derivation of the effective 
action for this quantum field theory is provided by writing down all terms in powers and 
gradients of $ which are invariant under all the symmetries of the Hamiltonian. In this 
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manner we obtain the action [3J H] 

= J d 2 rdr ((<9 T $) 2 + v 2 (V x <f>) 2 + s$ 2 + «($ 2 ) 2 ) (3) 

Here v is a spin-wave velocity, and s,u are parameters whose values adjusted to obtain 
Neel order in the ground state. In mean-field-theory, this happens for s < 0, where we 
have | (3?) | = (—s)/(2u) by minimization of the action S$>. A standard computation of 
the fluctuations about this saddle point shows that the low energy excitations are spin 
waves with two possible polarizations and an energy e which vanishes at small wavevectors 
k, e = vk. These spin waves correspond to local oscillations of $ about an orientation 
chosen by spontaneous breaking of the spin rotation symmetry in the Neel state, but which 
maintain low energy by fixing the magnitude |<fr|. The spinwaves also interact weakly with 
each other, and the form of these interactions can also be described by All effects of these 
interactions are completely captured by a single energy scale, p s , which is the 'spin stiffness', 
measuring the energy required to slowly twist the orientation of the Neel order across a large 
spatial region. At finite temperatures, the thermal fluctuations of the interacting spin-waves 
can have strong consequences. We will not describe these here (because they are purely 
consequences of classical thermal fluctuations), apart from noting jl] that all these thermal 
effects can be expressed universally as functions of the dimensionless ratio k^T j p s . 

For our future analysis, it is useful to have an alternative description of the low energy 
states above the Neel ordered state. For the Neel state, this alternative description is, in 
a sense, a purely mathematical exercise: it does not alter any of the low energy physical 
properties, and yields an identical low temperature theory for all observables when expressed 
in terms of fcgT '/ p s . The key step is to express the vector field $ in terms of a S = 1/2 
complex spinor field z a , where a =|j by 

$ = z* a a- al3 zp (4) 

where a are the 2x2 Pauli matrices. Note that this mapping from <fr to z a is redundant. 
We can make a spacetime-dependent change in the phase of the z a by the field 9(x, r) 

z a > e z a (5) 

and leave $ unchanged. All physical properties must therefore also be invariant under 
Eq. (pH), and so the quantum field theory for z a has a U(l) gauge invariance, much like that 
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found in quantum electrodynamics. The effective action for the z a therefore requires intro- 
duction of an 'emergent' U(l) gauge field (where \i = x, r is a three-component spacetime 
index). The field A^ is unrelated the electromagnetic field, but is an internal field which 
conveniently describes the couplings between the spin excitations of the antiferromagnet. As 
we have noted above, in the Neel state, expressing the spin-wave fluctuations in terms of z a 
and A^ is a matter of choice, and the above theory for the vector field $ can serve us equally 
well. The distinction between the two approaches appears when we move out of the Neel 
state across quantum critical points into other phases (as we will see later): in some of these 
phases the emergent A^ gauge field is no longer optional, but an essential characterization 
of the 'quantum order' of the phase. As we did for S$, we can write down the quantum 
field theory for z a and A^ by the constraints of symmetry and gauge invariance, which now 
yields 



S, = I d 2 rdr 



Kfy - iA^)z a \ 2 + s\z a \ 2 + u(\z a \ 2 ) 2 + ^(e^xduAx) 



2e 2 



(6) 



For brevity, we have now used a "relativistically" invariant notation, and scaled away the 
spin-wave velocity v ; the values of the couplings s, u are different from, but related to, those 
in S<f>. The Maxwell action for A^ is generated from short distance z a fluctuations, and it 
makes A^ a dynamical field; its coupling eo is unrelated to the electron charge. The action 
S z is a valid description of the Neel state for s < (the critical upper value of s will have 
fluctuation corrections away from 0), where the gauge theory enters a Higgs phase with 
(z a ) 7^ 0. This description of the Neel state as a Higgs phase has an analogy with the 
Weinberg-Salam theory of the weak interactions — in the latter case it is hypothesized that 
the condensation of a Higgs boson gives a mass to the W and Z gauge bosons, whereas here 
the condensation of z a quenches the A^ gauge boson. 



1. Triangular lattice 

There have been numerous recent studies [5] of the spin excitations of the insulator 
CS2CUCI4. Just as in La2Cu04, the dominant spin excitations are 5 = 1/2 spins on the 
Cu ions, but now they reside on the vertices of a triangular lattice, as sketched in Fig [TJd. 
Such an antiferromagnet is well described by the Hamiltonian H Q , with a nearest neighbor 
exchange J and i on the sites of the triangular lattice. From numerical studies of such spin 
systems [6], and also from observations [5] in CS2CUCI4, the ground state of H also has 
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broken spin-rotation symmetry, but the pattern of spin polarization is now quite different. 
We now replace Eq. ^ by 

(Sj) = Ni cos(K ■ rj) + N 2 sm(K ■ rj ) (7) 

where is the position of site i, and K = (47r/3a)(l, a/3) for the ordering pattern in Fig. [ljo 
on a triangular lattice of spacing a. The most important difference from Eq. ^ is that we 
now require two orthogonal vectors Ni^ (Ni • N2 = 0) to specify the degenerate manifold 
of ground states. As for the square lattice, we can write down an effective action for Ni^ 
constrained only by the symmetries of the Hamiltonian. Minimization of such an action 
shows that the ordered state has = N| fixed to a value determined by parameters 
in the Hamiltonian, but are otherwise arbitrary. Moving on to the analog of the spinor 
representation in Eq. Q, we now introduce another spinor w a , which parameterizes Ni^ by 

m 

Ni + lN 2 = (8) 

where e a p is the antisymmetric tensor. It can be checked that w a transforms as a S = 
1/2 spinor under spin rotations, and that under translations by a lattice vector y w a — > 
e~ lK ' y / 2 w a . Apart from these global symmetries, we also have the analog of the gauge 
invariance in Eq. (|5]). From the relationship of w a to the physical observables in Eq. ([8]), we 
now find a Z 2 gauge transformation 

w a -> T]W a (9) 
where f](r,r) = ±1. This Z 2 gauge invariance will play an important role in the discussion 



in Section IIP The low energy theory of the antiferromagnetically ordered state described 



by Eq. Q can now be obtained from the effective action for the Ni^ or the w a . We won't 



write it out explicitly here, deferring it also to Section IIP 



B. Coupled dimer antiferro magnet 

This spin model is illustrated in Fig. [2] We begin with the square lattice antiferromagnet 
in Figfjji, and weaken the bonds indicated by the dashed lines to the value Jjg. For g = 1, 
this model reduces to the square lattice model examined in Section II A| For g > 1, the 



model can be understood as a set of spin dimers, with the intra-dimer exchange interaction 
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FIG. 2: The coupled dimer antiferromagnet, described by the Hamiltonian Hq, with Jij = J on the 
full red lines, and Jij = J/g on the dashed red lines, (a) The large g ground state, with each ellipse 
representing a singlet valence bond — |l|))/\/2. (b) The 5 = 1 spin triplon excitation. The 
pair of parallel spins on the broken valence bond hops between dimers using the J/g couplings. 

J, and a weaker coupling between the dimers of J/g. A number of Cu compounds, such as 
TICUCI3 0,19] and BaCuSi206 [ID] , realize such coupled spin dimer models, although not in 
precisely the spatial pattern shown in Fig. [2} 

A simple argument shows that the large g ground state is qualitatively distinct from the 
g = 1 Neel state. For g = 00, the Hamiltonian separates into decoupled dimers, which have 
the obvious exact ground state 




where A represents the set of links with full lines in Fig. |2} This is the product of S = 
spin singlets on the dimers, and so preserves full rotational invariance, unlike the Neel state. 
There is a gap towards 5 = 1 excitations which are created by breaking the spin singlet 
valence bonds, as shown in Fig. [2Jd. This gap ensures that perturbation theory in 1/g is 
well-defined, and so the same structure is preserved for a finite range of 1/g. The most 
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significant change at non-zero 1/g is that the triplet excitations becomes mobile, realizing a 
S — 1 quasiparticle, sometimes called the triplon. Triplons have been clearly observed in a 
number of neutron scattering experiments on spin dimer [TT] and related [12] states. 

It is interesting to understand this state and its triplon excitations from the perspective 
of quantum field theory. By the same procedure as in Section II A we can argue that the 
theory in Eq. ^ applies also in the present phase. However, we now need s > to preserve 
spin rotation invariance. The excitations of S$> about the $ = saddle point consists of the 
3 polarizations of $ oscillations: we can identify these oscillations with a S — 1 quasiparticle, 
which clearly has all the characteristics of the triplon. 

So we now have the important observation that 5$ describes the g = 1 Neel state for 
s < s c (the critical value s c = in mean field theory), and the large g coupled dimer 
antiferromagnet for s > s c . This suggest that there is a quantum phase transition between 
these two states of H at a critical value g = g c which is also described by 5$. This is indeed 



the case, as we will discuss in Section III A 



C. Valence bond solids 



The state in Eq. ( |10[ ), and its triplon excitations, are also acceptable caricatures of a va- 
lence bond solid (VBS) phase. However, the crucial difference is that the VBS appears as the 
ground state of a Hamiltonian with full square lattice symmetry, and the space group symme- 
try is spontaneously broken, as shown in Fig. [3^. Thus Hq has all nearest-neighbor exchanges 
equal, and additional exchange couplings, with full square lattice symmetries, are required 
to destroy the Neel order and restore spin rotation invariance. Such VBS states, with spin 
rotation symmetry, but lattice symmetry spontaneously broken, have been observed recently 
in the organic compound [T3] EtMe 3 P[Pd(dmit) 2 ]2 and in [2] Zn a; Cu4_ ;c (OD) 6 Cl2. In the 
underdoped cuprates, bond-centered modulations in the local density of states observed 
in scanning tunneling microscopy experiments [15] have a direct interpretation in terms of 
predicted doped VBS states pH QH Q2] . 

We frame our discussion in the context of a specific perturbation, H 1 , studied recently 
by Sandvik [19] and others [201 EI]: 

Hi = -^(SrS r l) (S* - S, — |) , (11) 

(ijki) 
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FIG. 3: (a) Columnar VBS state of Hq + H\ with Jjj = J on all bonds. This state is the same as 
m Fig. |Ji, but the square lattice symmetry has been broken spontaneously. Rotations by multiples 
of 7r/2 about a lattice site yield the 4 degenerate states, (b) 4-fold degenerate plaquette VBS state 
with the rounded squares representing 5 = combination of 4 spins. 

where (ijkl) refers to sites on the plaquettes of the square lattice. The ground state of the 
Hamiltonian Hq + Hi will depend upon the ratio g = Q/J. For g = 0, we have the Neel 
ordered state of Section II A Recent numerical studies |19j [20], [21] have shown convincingly 



that VBS order is present for large g (VBS order had also been found earlier in a related 
'easy-plane' model [22J ) . We characterize the VBS state by introducing a complex field, 
whose expectation value measures the breaking of space group symmetry: 



(12) 



This definition satisfies the requirements of spin rotation invariance, and measures the lat- 
tice symmetry broken by the columnar VBS state because under a rotation about an even 
sublattice site by an angle = 0, n/2, n, 3tc/2 we have 



(13) 
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Thus \I/ is a convenient measure of the Z\ rotation symmetry of the square lattice, and (iff) = 
in any state (such as the Neel state) in which this Z4 symmetry is preserved. From these 
definitions, we see that (iff) 7^ in the VBS states; states with arg((\&)) = 0, ir/2, ir, 37r/2 
correspond to the four degenerate VBS states with 'columnar' dimers in Fig. |3^l, while states 
with arg((\l/)) = ir/i, 37r/4, 5n/A, 7tt /4 correspond to the plaquette VBS states in Fig. [3Jd. 

VBS order with 7^ was predicted in early theoretical work (23] on the destruction 
of Neel order on the square lattice. We now present the arguments, relying on the theory 
S z in Eq. (|6| as a representation of the excitations of the Neel state. We mentioned earlier 
that S z describes the Neel state for s < s c (with s c = in mean field theory). Let us 
assume that S z remains a valid description across a quantum phase transition at s = s c , 
and ask about the nature of the quantum state [21] for s > s c . A conventional fluctuation 
analysis about the z a = saddle point reveals two types of excitations: (i) the S = 1/2 z a 
quanta which realize bosonic quasiparticles known as spinons (see Fig[4|, and (ii) the 
'photon' which represents a spinless excitations we can identify with the resonance between 
different valence bond configurations [2S] . More precisely, the 'magnetic flux' determines 




FIG. 4: Caricature of a spin liquid state. The valence bonds are entangled between different 
pairings of the spins, only one of which is shown. Also shown are two unpaired 5 = 1/2 spinons, 
which can move independently in the spin liquid background. 

the relative phases of different possible pairings of the electrons into valence bonds. It 
would also appear that the s > s c phase of S z breaks no symmetries of the underlying spin 
Hamiltonian Hq + H±. However, there are well-known arguments [26] that this phase of S z 
breaks a 'hidden' or 'topological' symmetry. Let us define the topological current by 

J» = ^uxd u A x (14) 
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where e^ v \ is the totally antisymmetric tensor in 3 spacetime dimensions, so that Jt is the 
'magnetic' flux. The conservation law, <9^J M = 0, is a trivial consequence of this definition, 
and represents the conservation of flux in the s > s c phase of S c . However, by defining 

J, = d,(, (15) 

where ( is a free massless scalar field at low energies, we can also identify J M as the Noether 
current associated with the shift symmetry ( — > ( + constant. This shift symmetry is 
spontaneously broken in the s > s c phase of S z , and ( is the associated Goldstone boson. 
(Indeed, this Goldstone boson is, of course, just the photon which has only a single 
allowed polarization in 2+1 dimensions.) Nevertheless, these observations appear purely 
formal because they only involve a restatement of the obvious divergencelessness of the flux 



in Eq. (14), and the shift symmetry appears unobservable. 

The key point made in Ref. I2"3l was that the above shift symmetry is, in fact, physically 
observable, and is an enlargment of the rotation symmetry of the square lattice. They 
focused attention on the monopole operator, V, which changes the total A^ flux by 2tt 
(by the Dirac argument, because total flux is conserved, any such tunneling event must 
be accompanied by a thin flux tube carrying flux 2tt — this string is unobservable to the 



z a quanta which carry unit A^ charge). From Eq. (15), the flux density is measured by 
d t ( } and so its canonically conjugate variable under the Maxwell action in Eq. ^ is C/ e o> 
consequently, the operator which shifts flux by 2n is 

V = exp (i^f) ■ (16) 

In terms of the underlying spin model, this monopole is a tunneling event between semiclas- 
sically stable spin textures which differ in their 'skyrmion' number. By a careful examination 
of the phase factors [27J associated with such tunnelling events, it is possible to show that V 



transforms just like \1/ in Eq. (13) under a lattice rotation (and also under other square 
lattice space group operations). However, this is also the action of V under the shift sym- 
metry, and so we may identify the Z 4 rotation with a particular shift operation. We may 
also expect that there are additional allowed terms that we may add to S z which reduce 
the continuous shift symmetry to only a discrete Z4 symmetry, and this is indeed the case 



see Eq. (20)). With this reduction of a continuous to a discrete symmetry, we expect 



that the Goldstone boson (the photon) will acquire a small gap, as will become clearer in 
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Section III B These arguments imply that we may as well identify the two operators 



V ~ * 



(17) 



with a proportionality constant that depends upon microscopic details. The breaking of shift 
symmetry in the s > s c phase of S z now implies that we also have 7^ 0, and so reach 
the remarkable conclusion that VBS order is present in this phase. Thus it appears that 
S z describes a transition from a Neel to a VBS state with increasing s, and this transition 



will be discussed further in Section |IIIB| We will also present there a comparison of the 
numerical observation of VBS order in Hq + H\ with present theory. 



D. Z2 spin liquids 



The discussion in Section |II C| on the entanglement between valence bonds illustrated in 
Fig. [4] raises a tantalizing possibility. Is it possible that this resonance is such that we obtain 
a ground state of an antiferromagnet which preserves all symmetries of the Hamiltonian ? 
In other words, we restore the spin rotation invariance of the Neel state but do not break 
any symmetry of the lattice, obtaining a state called resonating valence bond (RVB) liquid. 
Such a state would be a Mott insulator in the strict sense, because it has an odd number 
of electrons per unit cell, and so cannot be smoothly connected to any band insulator. For 
phases like those described by S z , with gapped spinons and a U(l) photon, the answer is 
no: the monopole operator V, identified with the VBS order 1 i r , has long-range correlations 
in a state with a well-defined U(l) photon. 

The remainder of the next two sections will discuss a variety of approaches in which this 
monopole instability can be suppressed to obtain RVB states. 

The route towards obtaining a RVB state in two dimensions was presented in Refs. |29l [301 
[3T1 The idea is to use the Higgs mechanism to reduce the unbroken gauge invariance from 
U(l) to a discrete gauge group, Z 2 , and so reduce the strength of the gauge-flux fluctuations. 
To break U(l) down to Z 2 , we need a Higgs scalar, A, that carries U(l) charge 2, so that 
A — > e 2l9 A under the transformation in Eq. ([5j. Then a phase with (A) 7^ would break 
the U(l) symmetry, in the same manner that the superconducting order parameter breaks 
electromagnetic gauge invariance. However, a gauge transformation with 9 = n, while acting 
non-trivially on the z a , would leave A invariant: thus there is a residual Z 2 gauge invariance 
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which constrains the structure of the theory in the Higgs phase. 

What is the physical interpretation of the field A, and how does its presence characterize 
the resulting quantum state i.e. what are the features of this Z 2 RVB liquid which distin- 
guishes it from other quantum states ? This is most easily determined by writing down 
the effective action for A, constrained only by symmetry and gauge invariance, including its 
couplings to z a . In this manner, we expand the theory S z to S z + S\, with 

S A = J d 2 rdr [Kfy - 2iA fM )A a \ 2 + S\A a \ 2 + u\A a \ 4 - iA a e a pz* a d a z* p + c.c] (18) 

We have introduced multiple fields A , with a spatial index a, which is necessary to account 
for the space group symmetry of the underlying lattice — this is a peripheral complication 
which we will gloss over here. The crucial term is the last one coupling A a and z a : it 
indicates that A is a molecular state of a pair of spinons in a spin-singlet state; this pair 
state has a "p-wave" structure, as indicated by the spatial gradient d a . It is now useful to 
examine the mean field phase diagram of S z + <Sa, as a function of the two "masses" s and 
s. We have 2 possible condensates, and hence 4 possible phases [T8] : 

(i) s < 0, s > 0: This state has (z a ) ^ and (A) = 0. We may ignore the gapped A modes, 



and this is just the Neel state of Section II A 



(it) s > 0, s > 0: This state has (z a ) = and (A) = 0. Again, we may ignore the gapped A 
modes, and this is the VBS state of Section II C| 



(Hi) s < 0, s < 0: This state has (z a ) ^ and (A) ^ 0. Because of the z a condensate, this 
state breaks spin rotation invariance, and we determine the spin configuration by finding 



the lowest energy z a mode in the background of a non-zero (A) in Eq. (18), which is 

z a = (w a e^ r + e a(3 w}e-^ r ) /y/2, (19) 

with w a a constant spinor. Inserting Eq. (19) into Eq. Q we find that $ is space-dependent 
so that (Si) obeys Eq. Q with Ni j2 given by Eq. ^ and the wavevector K = (ir, n) + 2(A). 



Thus this state is the same as the coplanar spin-ordered state of Section II A 1 ! The Z 2 
gauge transformation in Eq. ^ is the same as the Z 2 9 = ix transformation we noted earlier 
in this subsection. 

(iv) s > 0, s < 0: This state has (z a ) = and (A) ^ 0. This is the Z 2 spin liquid (or Higgs) 
state we are after. Spin rotation invariance is preserved, and there is no VBS order because 
monopoles are suppressed by the A condensate, (\&) ~ (V) = 0. 
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We are now in a position to completely characterize this Z 2 spin liquid. One class of 
excitations of this state are the z a spinons which now carry a unit Z 2 "electric" charge. A 
second class are vortices in which the phase A winds by 2tt around a spatial point: these 
are analogous to Abrikosov vortices in a BCS superconductor, and so carry flux n. 
Further, because the A M flux can be changed by 27r by monopoles, only the flux modulo 2-7T 
is significant. Consequently there is only a single type of vortex, which is often called a vison. 
Two visons can annihilate each other, and so we can also assign a Z 2 "magnetic" charge to a 
vison. Finally, from the analog of the Abrikosov solution of a vortex, we can deduce that a 
Z 2 electric charge picks up a factor of (—1) when transported around a vison. This Z 2 x Z 2 
structure is the fundamental characterization of the quantum numbers characterizing this 
spin liquid, and its electric and magnetic excitations. An effective lattice gauge theory (an 
'odd' Z 2 gauge theory) of these excitations, and of the transition between this phase and 
the VBS state was described in Refs. [32] (a detailed derivation was published later [33]). We 
also note that Z 2 x Z 2 is the simplest of a large class of topological structures described by 
Bais and collaborators [34, 35J in their work on Higgs phases of discrete gauge theories. 

Since the work of Refs. EH] 130] EU the same Z 2 x Z 2 structure has been discovered in a 
variety of other models. The most well-known is the toric code of Kitaev [36], found in an 
exactly solvable Z 2 gauge theory, which has found applications as a topologically protected 
quantum memory [37]. An exactly solvable model of S — 1/2 spins in a Z 2 spin liquid phase 
was also constructed by Wen [38J. Senthil and Fisher [£§] have described Z 2 spin liquid 
states in easy-plane magnets and superconductors. Sondhi and Moessner [10] have studied 
quantum dimer models which also exhibit Z 2 spin liquid phases. Finally, Freedman et al. 
[4"T] have presented a classification of topological quantum phases using Chern-Simons gauge 
theory in which the Z 2 x Z 2 state is the simplest case. 



III. QUANTUM CRITICAL POINTS AND PHASES 

In Section [TTJ we saw several examples of insulating quantum magnets which were tuned 
between two distinct quantum phases by a coupling g. We will now turn our attention 
to the transition between these phases, and show that the critical points realize non-trivial 
quantum states, as was promised in Section[T| In Section HI B we will also describe interesting 



cases where such critical quantum states are found not just at isolated points in the phase 
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diagram, but extend into a finite critical quantum phase. 



A. Landau-Ginzburg- Wilson criticality 



We begin by considering the coupled dimer antiferromagnet of Section [II B| With increas- 
ing g, this model has a phase transition from the Neel state in Fig. [T^i to the dimer state in 
Fig, [2^i. We also noted in Section II B that this transition is described by the theory in 



Eq. ([3]) by tuning the coupling s. From numerical studies [2], we know that the transition 
in the dimer antiferromagnet occurs at g = 1.9107. . . , and the measured critical exponents 
are in excellent agreement with those expected from The quantum state at g = g c is de- 
scribed by a conformal field theory (CFT), a scale-invariant structure of correlations which 
describes a renormalization group fixed point (known as the Wilson-Fisher fixed point) of 
S$. The ground state breaks no symmetries, and there is a spectrum of gap less excita- 
tions whose structure is not completely understood: the spectrum is partly characterized 
by the scaling dimensions of various operators in the CFT, which can be computed either 
numerically or using a number of perturbative field-theoretic methods. 

We place this critical quantum state in the category of Landau- Ginzbug- Wilson (LGW) 
criticality, because the above understanding followed entirely from an identification of the 
broken symmetry in the Neel phase as characterized by the field and a subsequent 
derivation of the effective action for $ using only symmetry to constrain the couplings 
between the long-wavelength modes. This procedure was postulated by LGW for classical 
thermal phase transitions, and it works successfully for the present quantum transition. 

Striking experimental evidence for this picture has been found recently jl2] in the coupled- 
dimer antiferromagnet TICUCI3. One of the primary predictions of S^> is that the low-lying 
excitations evolve from 2 spin wave modes in the Neel state, to the 3 gapped triplon particles 
in the dimer state. The extra mode emerges as a longitudinal spin excitation near the 
quantum critical point, and this mode has been observed in neutron scattering experiments. 

B. Critical U(l) spin liquids 



Consider the phase transition in the model Hq + Hi, as a function of the ratio g = Q/J, 
in a model which preserves full square lattice symmetry. The crucial difference from the 
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model of Section III A is that there is only one 5 = 1/2 spin per unit cell, and so there is 



no trivial candidate for a S = ground state (like that in Eq. (10)) in the non-Neel phase. 
Indeed, this argument is sufficient to exclude S$ from being the appropriate quantum field 
theory for the transition, because has a trivial gapped ground state for s large. 

Instead, as we discussed in Sections II A and |II C[ the ground state of H + Hi evolved 
with increasing Q/J from the Neel state in Fig. [T^i to the VBS states in Fig. [3j We also 
argued that the low energy theory for the vicinity of this transition was given by the theory 
S z in Eq. @ of S = 1/2 spinons coupled to a U(l) gauge field A^. The Neel order parameter 
<fr was related to the spinons by Eq. Q, while the monopole operator for the A^ gauge field 
was identified with the VBS order parameter \1/ in Eq. (17). Notice that this low energy 
effective theory is not expressed directly in terms of either the Neel order $ or the VBS 
order This proposal therefore runs counter to the LGW procedure, which would lead 
directly for a theory of coupled <E», \I/ fields, like those presented in Ref. HH Instead, the low 
energy theory involves an emergent U(l) gauge field A^. 

Now we ask about the fate of this U(l) gauge field across the Neel- VBS transition. We 



noted in Section II A that A^ was quenched by the Higgs phenomenon in the Neel state, 



while in Section II C we pointed out that An was gapped in the VBS state by the reduction 



of the continuous shift symmetry of the field ( in Eq. (15) to a Z± lattice rotation symmetry. 
However, is it possible that A^ becomes gapless and critical in the transition from the Neel 
to the VBS state ? 

An affirmative answer to the last question would yield a most interesting U(l) spin liquid 
state [281 WQ . For A^ to remain gapless, it is required that the continuous shift symmetry 
in C, be preserved. The connection with the lattice rotation symmetry implies that 
the resulting U(l) spin liquid state actually has a continuous spatial rotation symmetry. 
In other words, the wavefunction describing the resonating valence bonds has an emergent 
circular symmetry, even though the spins on the valence bonds reside on the sites of a square 
lattice. Recent numerical studies by Sandvik [19] and others [21] on the model H + Hi have 
observed this remarkable phenomenon. They computed the histogram in their simulation 
for the VBS order \l/ not too far from the transition, and the results are shown in Fig. |5j 
The most striking feature is the almost perfect circular symmetry of the histogram; similar 
results have also been obtained in Ref. [2D and earlier for an easy-plane model [22]- This 
circular symmetry also means that the system has not made a choice between the VBS 
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FIG. 5: Histogram of the VBS order parameter defined in Eq. (12) in the numerical study [19J by 



Sandvik of Hq + H\. The circular symmetry is emergent and implies a Goldstone boson which is 
the emergent photon A^. 

orders in Figs. [3^l and [3]d, for which the distribution in Fig. [5] would be oriented along 
specific compass directions (N,E,S,W, and NE,SE,SW,NW respectively). These numerical 
results strongly support our claim that the theory S z is the proper low energy theory for the 
transition from the Neel to the VBS state, rather than the analog of the theory in Ref. H3J 
Let us restate the emergence of circular symmetry in more formal terms. The theory 
S z has a global SU(2) spin rotation symmetry, and also a global U(l) symmetry associated 
with the rotation in phase of the monopole V due to the shift symmetry in massless scalar 



( (noted below Eq. (15)). The underlying spin model H + H\ also has a global SU(2) spin 
rotation symmetry, but only a Z 4 lattice rotation symmetry. The simplest term we can add 
to S z which breaks its U(l) symmetry down to Z4 is 



S v = A J d 2 rdr [V 4 + c.c] 



(20) 



and a non-zero A gives a mass to the scalar (. The circular symmetry of Fig [5] implies that 
A scales under renormalization to small values jH [28] over the length scales studied by the 
simulation. 

A separate question is whether the Neel- VBS transition is second-order, with a critical 
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point where A scales all the way to zero. Numerical results are inconclusive at present 
[TT?| |2"IH |2"T| H3J HH] , with indications of both a second-order and a weakly first-order transition. 
In any case, it is already clear from the present results (see Fig [5J that there is a substantial 
length scale, much larger than the spin correlation length, over which monopoles described 



by Eq. (20) are suppressed, and the spin liquid theory of z a spinons and the gauge field 
applies. Indeed, the length scale at which ( acquires a mass, and the z a spinons are confined, 
is larger than all system sizes studied so far. 

A separate class of critical U(l) spin liquids appear in models in which the spinons are 
represented as fermions. These states of antiferromagnets were proposed by Rantner and 
Wen jl?] , building upon early work of Affleck and Marston [IE] , and a complete theory was 
presented by Hermele et al. [331 [50]. They are closely analgous to the U(l) spin liquid 
discussed above but with two important differences: (i) the bosonic z a spinons are replaced 
by massless Dirac fermions ip a , and (ii) they are postulated to exist not at isolated critical 
points between phases but as generic critical quantum phases. The latter difference requires 
that the CFT of these spin liquids has no relevant perturbations which are permitted by 
the symmetries of the underlying antiferromagnet. This has so far only been established in 
a 1/N expansion in models with SU(iV) spin symmetry. These spin liquids have also been 
proposed [51] as candidate ground states for the kagome lattice antiferromagnet observed 
recently in Zn rr Cu 4 _ a; (OII) 6 Cl2. 



IV. SUPERFLUIDS AND METALS 



Our discussion of quantum phases and critical points has so far been limited to examples 
drawn from quantum magnetism, in which the primary degrees of freedom are stationary 
unpaired electrons which carry spin S = 1/2. We have seen that a plethora of interesting 
quantum states are possible, many of which have found realizations in experimental systems. 
In this section we will widen the discussion to including systems which have also allow for 
motion of mass and/or charge. Thus in addition to insulators, we will find states which are 
metals or superfluids/superconductors. We will find that many phases and transitions can 
be understood by mappings and analogies from the magnetic systems. 

The simplest examples of a quantum phase transition involving motion of matter are 
found in systems of trapped ultracold atoms, which have been reviewed earlier [52]. Here 
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the degrees of freedom are neutral atoms, and so the motion involves mass, but not charge, 
transport. Bosonic atoms are placed in a optical lattice potential created by standing waves 
of laser light, and their ground state is observed as a function of the depth of the periodic 
potential. For weak potential, the ground state is a superffuid, while for a strong potential 
the ground state may be an insulator. 

A crucial parameter determining the nature of the insulator, and of the superfluid- 
insulator transition, is the filling fraction, /, of the optical lattice. This is the ratio of 
number of bosons to the number of unit cell of the optical lattice. 

We begin by considering the simplest case, / = 1, which is illustrated in Fig. [6] For a weak 



FIG. 6: Superfluid (top) and insulating states of bosons with repulsive interactions on a square 
lattice at / = 1. 

potential, the bosons move freely through the lattice, and form a Bose-Einstein condensate, 
leading to superfluidity. In the site representation shown in Fig. |6j the superfluid state is 
a quantum superposition of many particle configurations, only one of which is shown. The 
configurations involve large fluctuations in the number of particles on each site, and it is these 
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fluctuations that are responsible for easy motion of mass currents, and hence superfluidity. 
For a strong potential, the bosons remain trapped in the lattice wells, one each per unit 
cell, forming the simple 'Mott insulator' shown in Fig. [6} We can describe quantum phase 
transition between these states by using the LGW order-parameter approach, just as we 
did for the coupled-dimer antiferromagnet in Section |III A The order parameter now is the 
condensate wavefunction, which is represented by a complex field Indeed, this field has 
2 real, components but it is otherwise similar to the 3-component field $ considered in 

is then given by 



Section II A The LGW theory for the superfluid-insulator transition 
action in Eq. (J3l), with the replacement $ — > $. 

Now let us move to / = 1/2, which is illustrated in Fig. [7J The superfluid state is 



FIG. 7: Superfluid (top) and two possible insulating states of bosons on a square lattice at / = 1/2. 
The single bosons on the links in the first insulator resonate between the two sites on the ends of 
the link, and this insulator is similar to the VBS state in Fig. [3^,. Similarly the two bosons on the 
plaquette in the second insulator resonate around the plaquette, and this insulator is similar to the 
VBS state in Fig. Bb. 
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fairly similar to that at / = 1, with large number fluctuations on each site. However, the 
nature of the insulator is now very different. Indeed, in a single-band "boson Hubbard 
model" with purely on-site repulsive interactions between the particles, the ground state 
remains a superfluid even at infinitely large repulsion ('hard-core' bosons) : this is clear 
from Fig. [7j where it is evident that number fluctuations are possible at site even when 
double occupancy of the sites is prohibited. However, once off-site interactions are included, 
interesting insulating states become possible. The nature of these phases can be understood 
by a very useful mapping between hard-core bosons and the S = 1/2 antiferromagnets 



studied in Sections II A and II C ; the Hilbert spaces of the two models are seen to be identical 
by identifying an empty site with a down spin (say), and site with a boson with an up spin. 



Further, in this mapping, an easy-plane Neel state of Section II A maps onto the superfluid 
state of the bosons. It is then natural to ask for the fate of the VBS states in Fig. [3] under 
the mapping to the hard-core bosons. A simple analysis yields the insulating states shown in 
Fig. [7j The same reasoning also allows us to analyze the superfluid-insulator transition for 



f = 1/2 bosons, by connecting it to the Neel- VBS transition in Section III B In this manner, 
we conclude that the low energy quantum states near the superfluid-insulator transition are 
described by the U(l) gauge field theory S z in Eq. (|6]); the only change is that an additional 
'easy-plane' term is permitted with lowers the global symmetry from SU(2) to U(l). The z a 
quanta of this theory are now emergent 'fractionalized' bosons which carry boson number 
1/2 (in units in which the bosons illustrated in Fig. [7] carry boson number 1); a microscopic 
picture of these fractionalized bosons can be obtained by applying the above mapping to 
the spinon states in Fig. |4| 

The crucial result we have established above is that LGW criticality does not apply at 
/ = 1/2. By extensions of the above argument it can be shown that LGW criticality applies 
only at integer values of /. Superfluid-insulator transitions at fractional / are far more 
complicated. 

The above gauge theory for bosons at / = 1/2 has been extended to a wide variety of 
cases, including a large number of other filling fractions |54j . and models of paired electrons 
on the square lattice [551 EE] also at general filling fractions. The resulting field theories 
are more complicated that S z , but all involve elementary degrees of freedom which carry 
fractional charges coupled to emergent U(l) gauge fields. 

In the hole-doped cuprate superconductors, there are a variety of experimental indications 
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[57] of a stable insulator in some compounds at a hole density of x = 1/8. As we will discuss 
further in Section [VJ it is useful to consider a superfluid-insulator transition at x = 1/8, 
and build a generalized phase diagram for the cup rates from it. The low energy theory for 
such a transition is similar [56] to that for bosons (representing paired holes) at f — 1/16. 
We show the spatial structure of some of the large number of allowed insulating states that 
are possible near the onset of superfluidity in Fig. [8] these are analogs of the two distinct 
possible state for / = 1/2 in Fig. [7] It is our proposal [T7J [18], [51] that the bond-centered 
modulations in the local density of states observed by Kohsaka et al. [15] (see Fig. [10] below) 
can be described by theories similar to those used to obtain Fig. [8} 

Having described the superfluid-insulator transition, we briefly turn our attention to 



the analogs of the 'exotic' spin liquid states of magnets discussed earlier in Sections II D 



and |III[ The topological order associated with such states is not restricted to insulators 
with no broken symmetry, but can be extended not only into insulators with conventional 
broken symmetries, but also into superfluids [39]. Indeed, an exotic metallic state is also 
possible [58], called the fractionalized Fermi liquid (FFL). The simplest realization of a 
FFL appears in two band model (a Kondo lattice), in which one band forms a Z 2 spin 
liquid, while the other band is partially filled and so forms a Fermi surface of ordinary 
electron-like quasiparticles. One unusual property of a FFL is that the Fermi surface encloses 
a volume which does not count all the electrons: a density of exactly one electron per 
unit cell is missing, having been absorbed into the Z 2 spin liquid. There is, thus, a close 
connection between the quantum entanglement of a spin liquid and deviations from the 
Luttinger counting of the Fermi surface volume. The FFL also has interesting quantum 
critical points towards conventional metallic phases. These have been reviewed in Ref. 1591 and 
have potential applications to quantum phase transitions in the rare-earth 'heavy fermion' 
compounds. 

Another class of metallic and superconducting states, the 'algebraic charge liquids' [60] 
are obtained by adding mobile carriers to antiferromagnets described by <S Z . In this case, 
the charge carriers are spinless fermions which carry charges of the gauge field. These 
fermions can form Fermi surfaces or pair into a superconductor. These exotic metallic and 
superconducting states have a number of interesting properties which have been argued to 
describe the underdoped cuprates 
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FIG. 8: Structure of possible insulating states on the square lattice for paired electrons at a hole 
density of x = 1/8 (/ = 1/16) from Ref. O The shading is a measure of the charge density; the 
density is modified on bonds or plaquettes for configurations like those of the insulating states in 
Fig.0 



V. FINITE TEMPERATURE QUANTUM CRITIC ALITY 



We now turn to an important question posed in Section |Tj what is the influence of a zero 
temperature quantum critical point on the T > properties ? 

Let us address this question in the context superfluid-insulator transition of bosons on a 
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square lattice at / = 1 discussed in Section IV The phase diagram of this model is shown in 
Fig. [9] as a function of T and the tuning parameter g (the strength of the periodic potential). 
Note the region labeled 'quantum critical': here the physics is described by taking the critical 
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FIG. 9: Phase diagram of the superfluid-insulator transition in two dimensions. Similar crossover 
diagrams apply to most of the quantum phase transitions discussed in this paper. The quantum 
critical point is at g = g c , T = 0. The dashed lines are crossovers, while the full line is a phase 
transition at the Kosterlitz-Thouless temperature Tkt > 0. 



Hamiltonian at g — g c and placing it at a temperature T. The deviation of the coupling 
g from g c is a subdominant perturbation — upon examining the quantum fluctuations with 
decreasing energy scales, the excitations of the quantum critical state are first damped out 
by thermal effects, before the system is able to decide which side of the transition it resides 
on i.e. before it is sensitive to whether g is smaller or larger than g c . Thus we reach 
the seemingly paradoxical conclusion that the window of influence of the quantum critical 
point actually widens as temperature is raised. Of course, there is also an upper-bound in 
temperature, whose value depends upon all details of the particular experimental system, 
above which the quantum criticality is irrelevant because thermal effects are so strong that 
neighboring quantum degrees of freedom barely detect each other's presence. 

The 'transport' properties of systems in the quantum critical region are of particular in- 
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terest. For the superfluid-insulator transition, this means the study of the electrical conduc- 
tivity. More generally, we can consider the response functions associated with any globally 
conserved quantity, such as the total spin of the quantum magnetic systems of Section |TT} 
In traditional condensed matter physics, transport is described by identifying the low-lying 
excitations of the quantum ground state, and writing down 'transport equations' for the 
conserved charges carried by them. Often, these excitations have a particle-like nature, such 



as the 'triplon' particles of Section |II B[ in which case we would describe their collisions 
by a Boltzmann equation. In other cases, the low-lying excitations are waves, such as the 



spin- waves of the Neel states of Section II A and their transport is described by a non- 
linear wave equation (such as the Gross-Pitaevski equation). However, for many interesting 
quantum critical points, including the superfluid-insulator transition in two dimensions, nei- 
ther description is possible, because the excitations do not have a particle-like or wave-like 
character. 

Despite the absence of an intuitive description of the quantum critical dynamics, we 
can expect that the transport properties should have a universal character determined by 
the quantum field theory of the quantum critical point. In addition to describing sin- 
gle excitations, this field theory also determines the S-matrix of these excitations by the 
renormalization group fixed-point value of the couplings, and these should be sufficient to 
determine transport properties [61]. So the transport co-efficients, like the conductivity, 
should have universal values, analogous to the universality of critical exponents in classical 
critical phenomena. 

Let us now focus on critical points described by CFTs, and on the transport properties of 
CFTs at T > 0. CFTs in spatial dimension d = 1 are especially well understood and familiar, 
and so we consider e.g. the superfluid-insulator transition in one dimensional systems. A 
phase diagram similar to that in Fig. [9] applies: the quantum critical region of interest to 
us here is retained, but the critical temperature for the onset of superfluidity is suppressed 
to zero. We will compute the retarded correlation function of the charge density, x(k,&), 
where k is the wavevector, and u is frequency; the dynamic conductivity, ct{uj), is related to 
X by the Kubo formula, <j{uj) = limk^o(—iuj/k 2 )x{k,uj). It can be shown that all CFTs in 
d = 1 have the following universal x : 

4e 2 vk 2 4e 2 vK 

X(k,u) = -r-K—— - — , , <t(u) = -r-^—^ ! d = 1, all hu/k B T (21) 

n [v A k l — [io + irj) 2 ) h {—too) 



28 



where e and h are fundamental constants which normalize the response, K is a dimensionless 
universal number which characterizes the CFT, v is the velocity of 'light' in the CFT, and rj 
is a positive infinitesimal. This result holds both at T = and also at T > 0. Note that this 
is a response characteristic of freely propagating excitations at velocity v. Thus, despite the 
complexity of the excitations of the CFT, there is collisionless transport of charge in d = 1 
at T > 0. Of course, in reality, there are always corrections to the CFT Hamiltonian in any 
experimental system, and these will induce collisions — these issues are discussed in Ref. E21 
It was argued in Ref. EH that the behavior of CFTs in higher spatial dimension is qual- 
itatively different. Despite the absence of particle-like excitations of the critical ground 
state, the central characteristic of the transport is a crossover from collisionless to collision- 
dominated transport. At high frequencies or low temperatures, we can obtain a collisionless 



response analogous to Eq. (21), using similar arguments based upon conformal invariance, 

4 ^ /c^ 4 ^ 

x(k,u) = -^-K r== = , a(u) = ^-K ; d = 2,hu^>k B T, (22) 
h ^v 2 k 2 - (uj + zr]) 2 h 



where again K is a universal number. However, it was argued [61] that phase-randomizing 
collisions were intrinsically present in CFTs in d = 2 (unlike d = 1), and these would 
lead to relaxation of perturbations to local equilibrium and the consequent emergence of 
hydrodynamic behavior. So at low frequencies, we have instead an Einstein relation which 
determines the conductivity with [63] 

Dk 2 4e 2 
X (k,u) = 4e2 Xc^p— , = Ae\ c D = ; d = 2, hu < k B T, (23) 

where \ c is the compressibility and D is the charge diffusion constant of the CFT, and 
these obey \ c = QiksT/(hv ) 2 and D = &2hv 2 /(ksT) where 0^2 are universal numbers. A 
large number of papers in the literature, particularly those on critical points in quantum 



Hall systems, have used the collisionless method of Eq. (22) to compute the conductivity. 



However, the correct d.c. limit is given by Eq. (23); given the distinct physical interpretation 
of the two regimes, we expect that K ^ 0i©2> an d so a different conductivity in this limit. 
This has been shown in a resummed perturbation expansion for a number of CFTs [64] . 

Much to the surprise of condensed matter physicists [65] , remarkable progress has become 
possible recently on the above questions as a result of advances in string theory. The 
AdS/CFT correspondence connects the T > properties of a CFT to a dual theory of 
quantum gravity of a black hole in a negatively curved "anti-de Sitter" universe [66l EH 
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[70] : the temperature and entropy of the CFT determine the Hawking temperature 
and Bekenstein entropy of the black hole. This connection has been especially powerful 
for a particular CFT, leading to the first exact results for a quantum-critical point in 2+1 
dimensions. The theory is a Yang-Mills non-Abelian gauge theory with a SU(iV) gauge group 
(similar to quantum chromodynamics) and M = 8 supersymmetry. In 2+1 dimensions the 
Yang Mills coupling g flows under the renormalization group to a strong-coupling fixed point 
g = g c : thus such a theory is generically quantum critical at low energies, and in fact realizes 



a critical spin liquid similar to those considered in Section III B The CFT describing this 
fixed point is dual to a particular phase of M theory, and this connection allows computation, 
in the limit of large N, of the full x(A;, u) of a particular conserved density. These results 
contain [71] the collisionless to hydrodynamic crossover postulated above between exactly 



the forms for x(fc,u;) in Eq. (22) and (23), including the exact values [TTj, |T2] of the numbers 
K, Oi, 9 2 . A curious feature is that for this solvable theory, K = ©i©2, and this was 
traced to a special electromagnetic self-duality property in M theory, which is not present 
in other CFTs. It is worth noting that these are the first exact results for the conductivity 
of an interacting many body system (there are also conjectured exact results for the spin 
diffusivity of one-dimensional antiferromagnets with a spin gap [73], in a phase similar to 



Section II B whose dynamics was observed in Ref. fT2]) . 



A. Hydrodynamics 

The transport results for the superfluid-insulator transition presented so far have an 
important limitation: they are valid only at the specific particle densities of the insulators, 
where there are a rational number of particles per unit cell e.g. the insulators in Figs. |6j 
[7| and [8j While the particle density is indeed pinned at such commensurate values in the 
insulator at T = 0, there is no pinning of the density in the superfluid state, or anywhere 
in the phase diagram at T > 0. Even at small variations in density from the commensurate 
value in these states, we expect strong modifications of the transport results above. This is 



because Eqs. (22) and (23) relied crucially joT] on the lack of mixing between the conserved 
momentum current (in which the particle and hole excitations are transported in the same 
direction) and the electrical current (which is not conserved because of collisions between 
particles and holes moving in opposite directions). When particle-hole symmetry is broken 
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by a variation in the density, the mixing between these currents is non-zero, and a new 
analysis is necessary. 



As a particular motivation for such an analysis, we present in Fig. 10 a hypothetical phase 
diagram for the hole-doped cuprate superconductors as a function of a chemical potential, 
/i, which controls the hole density. We also have a tuning parameter, g, which controls the 




FIG. 10: Superfluid-insulator transition applied to the cuprate superconductors. The parameter 
g is as in Fig. [9j and that figure applies only exactly at hole density x = 1/8. The above figure 
shows how the quantum critical point (red circle) controls the phase diagram at general chemical 
potential fi and hole density x. We hypothesize that the density modulation in the insulator is 
similar to the bond-centered modulations in the underdoped superconductor observed by Kohsaka 
et al. PS], as indicated by the picture of their observations. Possible supersolid phases are not 
shown. 



proximity to an insulating state whose features have presumably been measured in Refs. [15] 
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and [571 We are interested in using the quantum critical point shown in Fig. 10 to describe 
transport in the analog of the quantum critical region of Fig. [9] when extended away to 
generic densities. Further, it is also interesting to add impurities and a magnetic field, B, 
to be able to describe various magnetotransport experiments |74j . 

Further motivation for the above analysis arises from potential applications to graphene 
|75j . The electronic excitations of graphene have a Dirac-like spectrum, and after including 
their Coulomb interactions, we obtain a system with many similarities to quantum critical 
states [7B]. At zero bias, we expect the conductivity of pure graphene to be described by 



Eqs. (22) and (23) (with logarithmic corrections from the marginally irrelevant Coulomb 
interactions [64J ) , but we are interested in the bias dependence of the conductivity, and also 
in the influence of impurity scattering. 

A general theory of quantum critical transport with all of the above perturbations has not 



been achieved. However, in the low frequency, collision-dominated regime of Eq. (23), it has 
recently been shown [77] that a nearly complete description can be obtained by using general 
hydrodynamic arguments based upon the positivity of entropy production and relaxation to 
local equilibrium. The results were verified in an explicit solution of the hydrodynamics of a 
particular CFT solvable via the AdS/CFT correspondence: after the addition of a chemical 
potential and a magnetic field, the CFT maps onto a black hole with electric and magnetic 
charges. 

The complete hydrodynamic analysis can be found in Ref. [77J along with a comparison 
to experimental measurements of the Nernst effect in the cuprate superconductors [74J; 
reasonable agreement is obtained, with only two free parameters. The analysis is intricate, 
but is mainly a straightforward adaption of the classic procedure outlined by Kadanoff and 
Martin [78J to the relativistic field theories which describe quantum critical points. We list 
the steps: (i) Identify the conserved quantities, which are the energy- momentum tensor, 
T^, and the particle number current, J M . (m) Obtain the real time equations of motion, 
which express the conservation laws: 

d v T^ = F» U J U , 9„J" = 0; (24) 

here F^ v represents the externally applied electric and magnetic fields which can change 
the net momentum or energy of the system, and we have not written a term describing 
momentum relaxation by impurities. (Hi) Identify the state variables which characterize 
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the local thermodynamic state — we choose these to be the density, p, the temperature T, 
and an average velocity u^ L . (iv) Express T^ v and J^ 1 in terms of the state variables and 
their spatial and temporal gradients; here we use the properties of the observables under a 
boost by the velocity and thermodynamic quantities like the energy density, s, and the 
pressure, P, which are determined from T and p by the equation of state of the CFT. We 
also introduce transport co-efficients associated with the gradient terms, (v) Express the 
equations of motion in terms of the state variables, and ensure that the entropy production 
rate is positive [79J. This is a key step which ensures relaxation to local equilibrium, and 
leads to important constraints on the transport co-efficients. In d = 2, it was found that 
situations with the velocity u M spacetime independent are characterized by only a single 
independent transport co-efficient [77J. This we choose to be the longitudinal conductivity 
at B — 0. (vi) Solve the initial value problem for the state variables using the linearized 
equations of motion, (mi) Finally, translate this solution to the linear response functions, 
as described in Ref. EHJ 

As an example of the above analysis, we highlight an important result-the modification 



of the conductivity in Eq. (23) by a non-zero density difference, p, from that of the Mott 
insulator at B — 0. We found that a non-zero p introduces an additional component of 
the current which cannot relax by collisions among the excitations, and which requires the 
presence of impurities to yield a finite conductivity: 

a (uj) = =-Q X Q 2 + y— — — -^7- -. — r. (25) 



Here the first term is as in Eq. (23), arising from the collisions between the particle and hole 
excitations of the underlying CFT, although the values of Q% and O2 can change as function 
of p. The additive Drude-like contribution arises from the excess particles or holes induced 
by the non-zero p, and l/r imp is the impurity scattering rate (which can also be computed 
from the CFT [7U ISO]). Note that in the 'non-relativistic limit', achieved when p is so 
large that we are well away from the critical point, the energy and pressure are dominated 
by the rest mass of the particles e + P ~ \p\mv 2 , and then the second term reduces to 
the conventional Drude form. The above theory can also be extended to include the Hall 
response in a magnetic field, B, [771 EE]; m particular we obtain the d.c. Hall conductivity, 
a xy oc p/B. This predicts a negative Hall resistance in the cuprates at hole densities smaller 
than that of the insulator at x — 1/8, as has been observed in recent experiments on the 
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cuprates [82] . 

VI. CONCLUSIONS 

We have described progress in understanding new types of entanglement in many elec- 
tron systems. Entanglement is a key ingredient for quantum computing, and the critical spin 
liquids described here are likely the states with most complex entanglement known. This 
progress has been made possible by advances on several distinct fronts. On the experimental 
side, new 'correlated electron' materials have been discovered, such as those described in 
Refs. [131 HH There have been advances in experimental techniques, such as scanning tunnel- 
ing microscopy in Ref. [151 an d neutron scattering in Refs. [121 And a new frontier has opened 
in the study of correlated phase of cold atoms in optical lattices |52j . On the theoretical side, 
large scale numerical studies of new types of quantum magnets has become possible by new 
algorithms, including the first studies of phases with no magnetic order in S = 1/2 square 
lattice antiferromagnets [191 EDI El] . An important frontier is the extension of such methods 
to a wider class of frustrated quantum systems, which have so far proved numerically in- 
tractable because of the 'sign' problem of quantum Monte Carlo. In field-theoretical studies 
of novel quantum phases, new connections have emerged to those studied in particle physics 
and string theory [5^1 \7T\ [77] . and this promises an exciting cross-fertilization between the 
two fields. 
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